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HYPOELLIPTIC MULTISCALE LANGEVIN DIFFUSIONS: LARGE 
DEVIATIONS, INVARIANT MEASURES AND SMALL MASS 

ASYMPTOTICS 

WENQING HU AND KONSTANTINOS SPILIOPOULOS 


Abstract. We consider a general class of non-gradient hypoelliptic Langevin diffu¬ 
sions and study two related questions. The first one is large deviations for hypoellip¬ 
tic multiscale diffusions. The second one is small mass asymptotics of the invariant 
measure corresponding to hypoelliptic Langevin operators and of related hypoelliptic 
Poisson equations. The invariant measure corresponding to the hypoelliptic problem 
and appropriate hypoelliptic Poisson equations enter the large deviations rate function 
due to the multiscale effects. Based on the small mass asymptotics we derive that the 
large deviations behavior of the multiscale hypoelliptic diffusion is consistent with the 
large deviations behavior of its overdamped counterpart. Additionally, we rigorously 
obtain an asymptotic expansion of the solution to the related density of the invariant 
measure and to hypoelliptic Poisson equations with respect to the mass parameter, 
characterizing the order of convergence. The proof of convergence of invariant mea¬ 
sures is of independent interest, as it involves an improvement of the hypocoercivity 
result for the kinetic Fokker-Planck equation. We do not restrict attention to gradient 
drifts and our proof provides explicit information on the dependence of the bounds of 
interest in terms of the mass parameter. 

Keywords: Large deviations, hypoelliptic multiscale diffusions, homogenization, 
hypocoercivity, non-gradient systems 
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1. Introduction 

The second order Langevin equation 

Tqt = f{qt) - Mt + cr{qt)Wt , go = O' e M” , go = P e M” , 

is one of the most classical equations in probability theory as well as in mathematical 
physics ( |18t fTTl 126]). It describes, under Newton’s law, the motion of a particle of 
mass T in a force field /(g), g G M"’, subject to random fluctuations and to a friction 
proportional to the velocity. Here Wt is the standard Wiener process (Brownian motion) 
in M”’, A > 0 is the friction coefficient. 

In this paper we are interested in the case where the force field /(g) has multiscale 
structure and the magnitude of the random fluctuations are small when allowing for 
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inhomogeneous friction coefficient. In particular, our starting object of interest is the 
second order hypoelliptic multiscale Langevin equation, 


( 1 ) 


rii! = 



+ c 



A {q!) qt 


dt + ^/ea 


qt 




where e, d <C 1 and 5 = (5(e) | 0 as e | 0. Here, X{q) > 0 is an inhomogeneous friction 
coefficient. Moreover, e represents the strength of the noise, whereas 6 is the parameter 
that separates the scales. 

It is well known that when r j, 0, the solution to ([T|) approximates that of a first 
order equation. In particular, if A is a constant, then in the overdamped case, i.e. when 
T is small, the motion can be approximated by the first order Langevin equation (see 
for example |12j i 


( 2 ) 
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The situation is much more complex in the case that the friction coefficient depends 
on the position too, see nails]. In particular, in the setting of ([T|), the motion of as 
r I 0 is approximated by 


( 3 ) 
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where a{q,r) = a{q,r)a'^{q,r). Clearly, when A(g) = A = constant, ([3|) reduces to ([2|). 

The first goal of this paper is to consider, in the homogenization regime where 
I —>■ oo as e, (5 j, 0, the large deviations behavior of the solution to ([I|) in such a way 
that, when the mass is small, it is consistent with the large deviations behavior of the 
solution to the overdamped counterpart Qj or equivalently ([2|). In particular we want 
to investigate the conditions under which the tail behavior of ([T]) and of ([3|) agree, at 
least in a limiting sense. 

It turns out that we get interesting non-trivial behavior when the mass r relates to 
e, (5 in a specific way that will be explained in the sequel. For this reason we shall write 
in place of r when we want to emphasize this dependence. We prove that if the mass 
of the particle r scales appropriately with the order of the fluctuations and in particular 
if it is of order (5^/e, i.e., if r = with m small but positive, then the large deviation 
behaviors of the overdamped and underdamped systems agree. The large deviations 
result for ([T]) is given in Theorem [2T] and the agreement in terms of the large deviations 
behavior of ([T|) and ([3|) is given in Theorem 12.51 

In order to derive the large deviations principle we follow the weak convergence 
approach, leiE]- This framework transforms the large deviations problem to convergence 
of a hypoelliptic stochastic control problem. Due to the hypoellipticity one needs certain 
a-priori bounds that establish compactness, see [T3|. We obtain an explicit form of the 

control (equivalently change of measure) that leads to the proof of large deviations upper 
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bound in the multiscale hypoelliptic case. Even though we do not address this issue in 
the current paper, we mention that the explicit information on the optimal control can 
be used for the construction of provably-efficient Monte Carlo schemes in the spirit of 
the constructions done in for the corresponding elliptic case. 

Under the parameterizations r = and when 5 <C e we derive the large devia¬ 
tions principle for {q^,e > 0}, where solves ([T|), see Theorem l2.ll The large deviations 
rate function is derived in closed form and it depends on m. The next natural question 
is to derive that as m | 0 the large deviations rate function converges to that of the 
large deviations principle for the overdamped case, i.e., for the solution to (l3|). However, 
to our surprise, we find that even in the case of constant diffusion the proof of such 
a convergence is highly involved. We prove such a convergence in the special case of 
diffusion coefficient a{q, r) = /3 > 0 (which is the parametrization of the 

fluctuation-dissipation theorem) and we include a discussion for the general variable 
diffusion coefficient case in Remark 12.61 This result supports the claim that the large 
deviations behavior of the multiscale second order Langevin diffusion ([T]) and of its first 
order counterpart dH) agree, see Theorem 12.51 

The second and related goal of this paper is to rigorously develop small mass 
asymptotics for the invariant measure, see Theorem 12.31 and for certain Poisson equa¬ 
tions, see Theorem 12.41 that appear in the rate function of the large deviation principle 
(see Theorem EH) due to the homogenization effects. Our proof of the convergence as 
m 4, 0 of the large deviation rate function requires a thorough analysis of the small mass 
asymptotic for the invariant measure of the fast motion corresponding to ([T])- In Par¬ 
ticular, since we will allow the drift term h{q, r) to be a general vector field rather than 
a gradient field, our proof of the convergence involves a non-trivial improvement of the 
hypocoercivity result for linear Fokker-Planck equation ([30l Section 1.7], see also [5]). 
If h{q, r) is not a gradient field, then certain operators that appear in the analysis are 
not anti-symmetric. This implies that extra terms appear that need to be appropriately 
handled. Then making use and extending the hypocoercivity results of m , we prove 
that the invariant measures corresponding to the m > 0 case, converges in to the 
invariant measure corresponding to the m = 0 problem. Here we make use of the ((•, •)) 
inner product introduced in [30] and we combine the different terms in such a way that 
the desired bounds follow. To accomplish this goal in the general non-gradient case, we 
use the structure of the hypoelliptic operator in an effective way. 

Using the convergence of the invariant measure and Poincare inequality, we also 
prove that the solutions to related Poisson equations (the so-called “cell problems”) 
that appear due to the homogenization effects of the drift b(q,r), also converge in the 
appropriate sense. In addition to that, the proof provides a rigorous justification of 
the corresponding multiscale expansion of the solutions of the corresponding equations 
in powers of ^/m. Related heuristic, i.e., without proof, asymptotic expansions can be 
also found in [25| . We would like to emphasize that our method of proof allows to obtain 


3 



upper bounds for the norms of interest with detailed dependence on the parameters of 
interest, such as the mass of the particle. 

Partial motivation for our work comes from chemical physics and biology, and in 
particular from the dynamical behavior of proteins such as their folding and binding 
kinetics. As it has been suggested long time ago (e.g., [laiM]) the potential surface 
of a protein might have a hierarchical structure with potential minima within poten¬ 
tial minima. As a consequence, the roughness of the energy landscapes that describe 
proteins has numerous effects on their kinetic properties as well as on their behavior at 
equilibrium. 

One of the first papers that used a simple model with two separated time scales to 
model diffusion in rough potentials is [55]. The situation usually investigated [niEaEi 
is based on the first order equation ([2|) even though the physical model and what is 
many times used in molecular simulations is the more complex second order Langevin 
equation that involves both position and velocity, see for example [20], and would also 
usually include more than two separated time scales. The usual choice of coefficients is 
X{q) = constant, b{q,q/6) = -^VQ{q/5) , c{q,q/S) = --^W{q) and a{q,q/S) = 
\/2/3A/, where kjs is the Boltzmann constant and T is the temperature, in such a way 
that the fluctuations-dissipations theorem holds. We remark here that our formulation 
for the large deviations result is general and includes the parametrization suggested 
by the fluctuation-dissipation theorem as a special case. Notice that the choice of the 
separable drift 

b{q,q/b) = -VQ{q/5), c{q,q/5) = -W{q) 

represents the motion of a massless particle in a rough potential £Q{q/5) + V{q). In 
particular, the model of interest in this case becomes 

(!) - «< > + . 

The questions of interest in [32l [9] are related to the effect of taking h | 0 with e 
small but hxed. This is almost the same to requiring that 5 goes to 0 much faster than 
£ does, which is the regime that we study in this paper. 

The related mathematical literature is quite rich. For the related hypocoercivity 
theory the reader is referred to m- For the case h = 1, the large deviations principle 
of the solutions to ([I]) and ([2]) as e 0 is being compared in [3]. For the case e = 1, 
periodic homogenization for a special case of ([T]) (in particular when c{q, r) = 0 and 
b{q,'£) = b{r)) has been addressed in [14]. Also, when e = 1 random homogenization for 
([T]) when c{q,r) = 0 and the special case of gradient drift b{q,r) = —VQ(r) has been 
addressed in [2l[23]. More is known about the over damped case ([2]), see [3 dzi [m |28] 
where homogenization and large deviation results for the solution to equations of the 
form ([2]) are obtained under different relations between e and 5, in both periodic and 
random environments. 
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The rest of the paper is structured as follows. In Section [2] we formulate the 
problem, our assumptions and the main results of this paper in detail. In Sections [3][3] 
we prove the large deviations principle for the hypoelliptic problem. In Sections [SHT] and 
in the Appendix we exploit the small mass asymptotics. 

In particular, using the weak convergence approach we turn the large deviations 
principle into a law of large numbers for a stochastic control problem. Section [3] proves 
the convergence of the controlled stochastic equation and Section [J] proves the conver¬ 
gence of the cost functional, which is the Laplace principle. In Section [5] we prove the 
small mass limit of the rate function in the diffusion a{q,r) = ^J2|3\{q)I case, using 
the convergence of the invariant measures as m —0 (Section [6]) and of the related “cell 
problems” that are auxiliary Poisson equations that appear in the rate functions due 
to homogenization effects (Section [7]). We emphasize that Section [6] is of independent 
interest as it is an extension of the hypo-coercivity result for the linear kinetic Fokker- 
Planck equation [30l Section 1.7], since we do not restrict our attention to drifts that 
are of gradient form. The method of proofs also yields explicit decay rates of the norms 
of interest with regards to parameters of interest such as the mass of the particle. Most 
of the proofs to technical lemmas are deferred to the Appendix. 


2. Problem formulation, assumptions and main results 

In this section, we formulate more precisely the problem that we are studying in 
this paper, we state our main assumptions and our main results. In preparation for 
stating the main results, we recall the concept of a Laplace principle. 


Definition 2.1. Let > 0} be a family of random variables taking values on a 

Polish space S and let I he a rate function on S. We say that {q^,£ > 0} satisfies 
the Laplace principle with rate function I if for every bounded and continuous function 
h : 5 —> M 


lim — elnE 
etO 


exp 



inf [L{x) + h{x)] . 


If the rate function has compact level sets, then the Laplace principle is equivalent to 
the corresponding large deviations principle with the same rate function (see Theorems 
2.2.1 and 2.2.3 in P). Hence, instead of proving a large deviations principle for we 
prove a Laplace principle for {q^}- 

Our main regularity assumption in regards to the coefficients of ([T]) is given by 
Condition 12.11 


Condition 2.1. The functions b{q,r),c{q,r), a{q,r) are 

(i) periodic with period 1 in the second variable in each direction, and 

(ii) C'^(M'^) in r and C'^(R'^) in q with all partial derivatives continuous and globally 
bounded in q and r. 
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The diffusion matrix a{q,r) = a{q,r)a'^{q,r) is uniformly non-degenerate. There exist 
constants 0 < A < A such that for every q £ A < \{q) < A. Moreover, the function 
X{q) is in with hounded partial derivatives. 

Using the parametrization r = the system being considered is 


( 5 ) 


... 

m —m = 
e 




+ 1 9(1 ^ 


Setting pf = y/m^ql we obtain the following system of equations which we also 
supplement with initial conditions 

(6) « = 1 


m 5 


Pt = 


1 1 


m S 


. Qt 




. Qt 




m 


= go G , po = Po G 


Condition EH guarantees that dS]) and have a unique strong solution; this is a 
classical result, see for example m or Theorem 5.2.1 of [23]. The infinitesimal generator 
for the (g,p) process satisfying (|6|) is given by 


C = 


q/^) ■ '^p + ^^(q^ q/^) ■ ^p 


+ 


1 e 

m 5“^ 


-X{q)p . Vp + -a(g, q/6) : V: 


where we recall that a{q,r) = a{q,r)a'^{q,r). 

We can assume that po is a random variable, as 

^ UmaxPol 


62' 


ong as it is independent of the 
< oo (see Appendix]^, where 


driving Wiener process Wt and as long as E 
we have defined Umax = max sup \ ai^j{q,r) 

Sometimes, we may write Xf = {qf,Pt)- Let | • | be the Euclidean norm in M'^. We 
introduce the control set 

^ = |u = {u, £ : 0 < s < T} progressively ^'.-measurable and E^ < oo| , 

The result in [3] gives the following representation 

= inf Eg i / \us\‘^dsh{ql) 
ugA [27o J 
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Here the process gf is the (/“Component of the hypoelliptic controlle(i diffusion 
process Xf = {qf,pi): 


(7) 


1 e 


m 6 


Qt = 


1 1 

Pt = 


iPt 


m S 




V‘’ 

ro Pt ' jr ! 

m 0 ^ d \/m 


+ ^ ^ J ' Wt 


m 


^0 = e K"*, ^ = Po e K"* 


Let u% ^ A and solves dT]) with u’^ in place of u. Let the control space be 
Z = the fast variable space be 3^ = x T*^. We see that the fast variable is actually 

/_ q^\ 

{Psi~^j- us define the operator 

r) = \p ■ Vr<^>ip, r) + b{q, r) ■ Vp^>(p, r)] 

-Xiq)p ■ Vp<^{p, r) + ^a{q, r) : Vl^{p, r) . 

For each fixed q, the operator defines a hypoelliptic diffusion process on (p, r) e 3^ = 
X T'^. Let p{dpdr\q) be the unique invariant measure for this process. Notice that 
is effectively the operator corresponding to the fast motion. The following centering 
condition is essential for the validity of the results. 



Condition 2.2. ITe assume that for every q 

/ b{q,r)ii{dpdr\q) =0. 

Jy 

Let us consider the preliminary cell problem 

(8) £™$(p,r) =- ]=p , [ ^{p,r)p{drdp\q) = 0 . 

Vm Jy 

It is clear that the solution to ([8]) depends also on q, but we sometimes suppress 
this in the notation for convenience. By the work of [13], we know that under Condition 
[221 the PDF dH]) has a unique, smooth solution that does not grow too fast at infinity, 
see Appendix [XI for more details. Note that the function is actually a vector valued 
function $(p,r) = ($i(p, r),..., ^>rf(p, r)). 

Then our first main result reads as follows. 
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Theorem 2.1. Let {{q‘^,p^),e > 0} be the unique solution to Under Conditions \2.1 
and \2.2l {q^,e > 0} satisfies the large deviations principle with rate function 



- rmifis)) Qm {fis){4>s - r^{(f)s))ds iffi€ AC{[0,Ty,'R^),(f)o = qo 

otherwise . 


rm{q) = 


1 


m 


A 


Vp<I>(p, r)c{q, r)p{dpdr\q) 


Qmiq) — 


m 


'y 


Vp^{p, r)a{q, r)(Vp^>(p, r)f ij.{dpdr\q) . 


To support the claim that the particular parametrization is consistent with the large 
deviations principle of the overdamped case ([3]), we need to prove that lim Sm{4>) = 

m^O 

So{(j)), where So{(j)) is the rate function associated to (l3|). To that end, we recall the 
corresponding large deviations result from [7]. 

Let po{dr\q) be the unique invariant measure corresponding to the operator 


= 


b{q, r) ■Vr+ r) : 


2X{q) 


equipped with periodic boundary conditions in r (q is being treated as a parameter 
here). By Theorem 12.31 Condition 12.21 implies that the following centering condition for 
the drift term b: 

b{q,r)po{dr\q) = 0, 


ly 


where T = denotes the d-dimensional torus. Under this centering condition, the cell 
problem 


(9) 


^qXiiq, r) = ’ J- r)po{dr\q) = 0 ,£ 


= l,2,...,d 


has a unique bounded and sufficiently smooth solution x = (xit Xd)- After these 
definitions we recall the result from [7] that will be of use to us. 

Theorem 2.2 (Theorem 5.3 in [7]). Let {q^,£ > 0} be the unique solution to (0). 
Under Conditions \2.1\ and \2.A {q^,£ > 0} satisfies a large deviations principle with rate 
function 

- rQ{<j)s)Y'Qfi^{fis){fis - ro{(ps))ds if & .4C([0, T]; M'^), c/iq = qo 

otherwise. 



and 










In order now to show that lim Sm{4>) = 'S'o(<^)j we need to study the limiting 
begavior of ^{dpdr\q) and of Vp<h(p, r) as m —)• 0. For this purpose, let us assume 
that (7{q,r) = ^J2|3\{q)I, (3 > 0, i.e., we assume that the noise is such that we are in 
fluctuation-dissipation balance. In this case, for a function / £ C'^{y), we have 

= ^^Af{p,r) + -^Bf{p,r) , 

^ m y/m 

where Af = —p ■ Vpf + (3A.pf and Bf{p, r) = p ■ V^f + b{q, r) ■ Vpf ■ Likewise, we have 

= -;^h{q, r) ■ V^/(r) /3A^/(r) . 

Kq) 

We denote by p{dpdr\q) = p'^{p,r\q)dpdr the invariant measure corresponding to 
the operator Also, let us write po{dr\q) = pQ{r\q)dr for the invariant measure 

corresponding to the operator 

Let us also define 7r{dp) = p^^{p)dp to be the invariant measure on for the 
Ornstein-Uhlenbeck process with generator A. With this notation, let us write p^{p, r) = 
p^{p,r)p^{p,r), where p^{p,r) = p^^{p)po{r), suppressing the dependence on q. 

Then, in Sections E] and [3 respectively we prove the following Theorems which 
constitute the second main result of our paper. 


Theorem 2.3. Let Condition \2.1\ hold and assume that a{q, r) = y/2p\{q)I, (5 > 0. 
Then, for every 5 G we have 


hm ||;5™(p,r)-I 

m—^O 


L^{y\P°) 


= 0 . 


Theorem 2.4. Let Conditions \2. l\ and \2.2\ hold and assume that a{q,r) = y/2j3X{q)I, /3 > 
0. Then, for every g G we have 


lim 

m—>0 



1 

Kd) 


(/ + v,x) 


= 0 . 

L^[y\p°) 


Using then Theorems l2.3l and l2.4l we prove in SectionEJthat the rate function Sm{4>) 
converges So{fi), as m 0. 


Theorem 2.5. Let Conditions \2.1\ and \2.2\ hold and assume that a{q,r) = y/2ftX{q)I, fd > 
0. Then, we have 

lim Smif)) = So{(l)) . 
m^O 

Remark 2.6. We believe that Theorems I2.,‘?l and \2.4\ cind as a consequence Theorem 

I A 51 are true under more general variable diffusion coefficients as long as Condition \2.1\ 

holds. When, the diffusion coefficient a is not a multiple of the identity matrix, then 

the operator A is not the classical Ornstein-Uhlenbeck that has the Gaussian measure 
_\p\f 

p^^{p)dp e 2/3 dp ag itg invariant measure. Some of our technical lemmas use this 
explicit structure in order to derive the necessary estimates. However, since the spirit of 

the proof does not rely on this structure, we believe that this is only a technical problem. 
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3. Law of large numbers 


In this section we study the limiting behavior of the solution to the control problem 
Q . It turns out that we need to consider the solution to ([71) together with an appropriate 
occupation measure and then consider the limit of the pair. Let us be more specific now. 

Let u% ^ A and solves dTj) with in place of u. Let the control space be 
and the fast variable space be 3^ = x T'^. We see that the fast variable is actually 

P’1, . Let A Z, Bi X B 2 Q y and L C [0, T]. Let A = A(e) > 0 be a separation 

of scales parameter. We introduce the occupation measure 

f-t+A 


(10) X X ^2 X r) = ^ 

Let us dehne the function 

(11) j{q,{p,r),z) = 


-[ 


1a{uI)1bApI)^B2 


— mod 1 I ds 


dt . 


1 


m 


[c{q, r) + a{q, r)z] ■ Vp$(p, r) . 


Definition 13.11 captures the notion of a viable pair as introduced in [7] which char¬ 
acterizes the required law of large numbers. 


Definition 3.1. A pair {'tp, P) G C([0, T]; M“) xV{Z x 3^ x [0, T]) will be called viable with 
respect to ( 7 , £™) or simply viable if there is no confusion, if the following are satisfied. 
The function fit is absolutely continuous, P is square integrable in the sense that 


( 12 ) 

and 


L 


Zxyx[0,T] 


\z\'^P{dz,dpdr,ds) < 00 , 


(i) 

(13) fi^ = q^-y / 'y{fi^,{p,r),z)P{dz,dpdr,ds) ; 

Jo Jzxy 

(ii) For any g{p,r) G V^CAfi), 

(14) f [ Cfp^g{p,r)P{dz,dpdr,ds) =Q \ 

Jo Jzxy 

(hi) 


(15) P{Z xy x[0,t]) =t. 

We write {fi,P) G y{^,cTi- 

Theorem 3.1. Consider any family {u^,e > 0} of controls in A satisfying 

fT 

supE / lu’ifidt < 00 . 

£>0 Jo 

Let Conditions \2.1\ and \2.A\ he satisfied. Then the family > 0} is tight. 

Hence, given any subsequence of {{ql,P^’‘^),£ > 0}, there exists a subsequence that 
converges in distribution with limit (g,,P). With probability 1, the accumulation point 

( 7.5 P) 'Is a viable pair with respect to {'y,Cq): {q„P) G 
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Proof of Theorem \3.R Part 1. 
apply Ito’s formula and get 


[Tightness]. 


For a smooth function g 



we can 


g 





+ 




+ <7 






ds 


(16) 



dWs . 


Let us apply Ito’s formula to $ 
tation formula for qf as follows: 



in ([8]) and we use (fT6]l to get a represen- 



Using this representation formula, Condition 12.II and Theorem 3.3 of [H] (see also 
Appendix El), we can then establish that for every g > 0 


limlimsupP sup — V 

PtO etO \ti-t2\<p,0<ti<t2<l 


= 0 . 


This implies the tightness of the family {ql}- Tightness of the occupation measures 
> 0} follows from the bound 

rT 1 rt+A 

(17) sup E [^'(P'^’^)] = sup E / — \u^{s)\'^dsdt < oo. 

ee(0,l] £6(0,1] Jo ^ Jt 


for the tightness function g{r) = /^xy’xfoT] \^\'^'>’{dz,dpdr,dt), r G V{Z x y x [0,T]), 
see Theorem A. 19 in [6]. Notice that the last inequality in ()17p follows by the uniform 
bound on the family of controls {n^,e > 0}. 

Hence, the family {{q‘^, P*^’^), e > 0} is tight. Due to tightness, for any subsequence 
of e > 0 there exists subsubsequence that converges, in distribution, to some limit {q, P) 
such that 

(r,p^’^)^(g,p). 


Next, we prove that any accumulation point will be a viable pair according to Definition 

EH 

Part 2. [Proof of (|12p ]. By Fatou’s Lemma we have 


E [ 

Jzxyxio,T] 


z\‘^P{dz,dpdr,dt) < oo , 


which then implies that | 2 :|^P(dz, dpdr, dt) < oo w.p.l. 
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Part 3. [Proof of ([T^ ]. Consider a test function / = f{q) on Let '^{p,r) = 
^{p,r) ■ ^qf{q) which satisfies the cell problem 

/:^'L(p,r) = -^p.V,/(g) . 

Making use of ([8]) and (fT6l) we get 


m 


'I'is.f)-* p;.|i = - 


- f— 

<5^ Jo \M 


• ^qf{ts)ds 


+ 


1 /■* 1 


■L 


0 








(18) 


+ 


f\ 

6 \/m 


r 


Vp'h MIP. 


r 


Let us now choose S,t > 0 such that S < S + t < T. We have 

fiqs+r) - fi'fs) = [ Vg/(gf)dt . 

Js wrnd 

Combining the latter expression with (I18p we get 
fits+r) - fi'fs) - ■ ^qfiqt)dt 


= -6 ^ 


-fA 




m 


+ V^^I Vp^{pl,^]-a{q%^]dWs. 


-e r 


Due to the a-priori bounds from Appendix|^the right hand side of the last display 
goes to zero in L^, which means that 


fC^s+r) - f{qs) - 7 {qu {pu y) , • ^qf{qt)dt 

as e 4, 0 in means square sense. By Condition 12.11 Lemma 3.2 of [7] guarantees that 

7 (pt^ ■ ^qf{tt)dt 

-L 


2 x:yx[s,s-hr] 


and 


L 


.Ex:yx[ 5 ,S-l-r] 


7 ( 9 ?, (p, r), z) ■ Vgf{ql)F^’^{dz, dpdr, dt) 


7 ( 9 ?, (P, r), z) • Vq/(gf)P^’^(dz, dpdr, dt) 


■L 


Jzxyx[s,s+T] 
as e 4 , 0. Therefore, by dehning 

Te.A 


7(jo (P, r), z) ■ Vqf{q^)F{dz, dpdr, dt) 


0 


Af 


f(q)= [ j{q,{p,r),z)Vf{q)Ff‘^{dz,dpdr) , 
Jzxy 
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where 


Fl’‘^{dz,dpdr) = 


rt+A 


A 


Uz{ul)ldp {pD Idr -f mod 1 ds , 


we get that, as e J, 0, 


(19) 

and, in probability, 

rS+T 


E 


f{fs+r) - fifs) - [ 
Js 




( 20 ) 


Al 


f{fs)ds- 


zxyx[s,s+T] 


liQs, {p, r), z)Vf{qs)P{dz, dpdr, ds) 0. 


Relations (USD and (I20p imply that the pair (g, P) solves the martingale problem 
associated with (fT^ . which then proves that (fT^ holds. 


Part 4. [Proof of (1141) ]. For functions / € 6*^(3^), let us introduce the auxiliary 
operator 


1 1 


■K,qfiP^ 0 = ^) > 


and define the J^-martingale 


A/f = /(S.f )-/lr„,f ) - / Al..T.f[T< 


s' 




ds 


y/g 1 

6 Am 


f 


As 




Let us furthermore set Gf, zfip,r) = —= [c{q,r) + a(q,r)z] ■ Vpf{p,r) and define 

Wm 


g{e) = Then, we have that 
(21) g{e)Mf - g{e) 


+ 5(e) 


fM 


A 
5 

s+A 




■^U.rJ ( Po. % I 


I 


ds - I ■Ale,qff ( P. 


[P%^]- ( P% ^ 


dp 


^6 A 

ds 


ds 


Gqi,zf (p, r) P^’^(dz, dpdr, ds) 


g \Jzxyx[o,t] 

rt 2 fs+A 

Jo a / 




ffi" 


ds 


L 


TS/ (p, r) dpdr, dt). 


Zxyx[o,t] 
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Let us now analyze the different terms in (I2ip . We start by observing that E < 

which then implies that g{e)M^ 0 in probability, as e 4- 0- Moreover, bounded¬ 
ness of / implies that g{e) f converges to zero uniformly. Hence, 

the left hand side of (f^ converges to zero in probability as e | 0. 

Let us next study the right hand side of (1211) . We have the following 

(i) Conditions l2.ll the uniform bound on the controls and tightness of {g^, e > 0}, 
imply that the first and the third term in the right hand side of (1211) converge 
to zero in probability as (5/e I 0. 

(ii) The second term on the right hand side of (I2ip also converges to zero in prob¬ 
ability, by the fact that (5/e | 0 and uniform integrability of 

Thus, by combining the behavior of the different terms on the left and on the right 
hand side of (|2ip . we obtain that we should necessarily have that 

/ C’2f(j>,r)P'^’^{dz,dpdr,dt)^0, in probability. 

Jzxyx[o,T] ‘ 

which by continuity in t G [0,T] gives (1141) . 

Part 5. [Proof of (fTSl) .] Finally P{Z x y x [0, t]) = t follows from the fact that 
analogous property holds at the prelimit level, P{Z x y x {t}) = 0 and the continuity 
of t —)■ P{Z X T X [0, t]) and (fT5|) follows. □ 

4. Laplace principle 

The main result of this section is the following Laplace principle. During the proof 
of Theorem 14.11 we also establish the alternative representation of Theorem 12.11 

Theorem 4.1. Let {ql,e > 0} be the unique strong solution to (1). Assume Conditions 
\2.1\ and WM . Define 


( 22 ) 


Sm{4>) = , inf 

((i>,P)eV(^_£m) 


1 


\z\^P{dz, dpdr, dt) 


}Zxyx[0,T] 

with the convention that the infimum over the empty set is oo. Then for every hounded 
and continuous function h mapping C([0,T];M'^) into M we have 

KqIY 


lim — elnE, 
elO 


qo 


exp 


inf [Smif)) + h{(p)] . 
<7ieC([0,T];Kd) 


Moreover, for each s < oo, the set 

‘Ps = WeC{[0,T]-R'^):Sm{(^)<s} 

is a compact subset o/C([0,T];M'^). 

In other words, {g^,e > 0} satisfies the Laplace principle with rate function S{* 


Proof of Theorem \4-l\ The proof of this theorem borrows some of the arguments of the 
related proof of the LDP for the elliptic overdamped case of Theorem 2.6 in [7]. We 
present here the main arguments, emphasizing the differences. 
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Part 1. [Laplace principle lower bound]. Theorem l3.1l and Patou’s lemma, guarantee 
the validity of the following chain of inequalities. 


liminf ( —elnE 


exp 


> liminf ( E 




— S 


> liminf ( E 
etO 


= lim inf E 
etO \ 


1 /■*+^ 


1 r 1 r 

2 Jo a/ \ul\^ dsdt + hiq^^) 
zJP^'^{dz,dpdr,dt) + h{(f) 


IZxyx[0,T] 


> E 


'Zxyx[0,T] 

f 1 


zf P{dz, dpdr, dt) + h{q) 

zj^ P(dz, dpdr, dt) + h{(j)) 


— 

(?i,P)eV(;^_£m) Jzxyx[0,T] 

= inf [Sm{(f>) + h{(f))] . 

</>eC([o,r];K‘i) 

Hence, the lower bound has been established. 

Part 2. [Laplace principle upper bound and alternative representation]. We first 
observe that one can write (I22|l in terms of a local rate function 

r-T 

' ' is . 


Sm{ct>)= ! U{cj>s,^s)ds 
Jo 


Here we set 


L’'(x, v) = inf 


Jzxy 2 


}-\zfP{dz,dpdr) , 


where 


= < 


PeViZxy): [ C^fip, r)P{dz, dpdr) = 0, V/ G CL(3^) , 

Jzxy 

/ | 2 |^P((iz, dpdr) < oo and = / 'y{q,{p,r), z)P{dz,dpdr) 

Jzxy Jzxy 

We can decompose the measure P € V{Z x y) into the form 
P{dz,dpdr) = p{dz\p,r)p{dpdr\q) , 


where /r is a probability measure on y and p is a stochastic kernel on Z given y. This 
is referred to as the “relaxed” formulation because the control is characterized as a 
distribution on Z (given q and (p, r)) rather than as an element of Z. We now have, for 
every / G Cioc(3^) for every q G M'’*, that 



f{p,r)p{dpdr) = 0 . 


Here we have used the independence of T™ on the control variable z to eliminate 
the stochastic kernel p. Thus p(dpdr) is the unique corresponding to the operator C^, 
written as p{dpdr\q). 
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Since the cost is convex in z and 7 is affine in z, the relaxed control formulation is 
equivalent to the following ordinary control formulation of the local rate function 

1 


L°{q,v)= inf \v{p,r)\^^i{dpdr) , 


where 


■^q,u = ^ V 


{•) -.y V{y) , {v,p) satisfy [ f{p,r)p{dpdr) = 0,V/ e Cfoc(3^) > 

Jy 


/ |u(p, r)|^/i((ipdr) < 00 and = / 'y{q, {p,r),v{p,r))p{dpdr) 

Jy Jy 

One can show as in [71 Section 5] that U{q,v) = L°{q,ii). Let us recall now the 
dehnitions of rm{q) and Qm{q) from Theorem 12.11 For any v G we can write 


^ = i l{qAP^r),v{p,r))p{dpdr\q) 

1 ■ 

Jy 


= [c{q, r) + a{q, r)v{p, r)] ■ Vp$(p, r)p{dpdr\q) 

= rm{q)+ I -^Vp<^{p,r)a{q,r){v{p,r))'^p{dpdr\q) . 

Jy y/m 

Then, u — r^iq) can be treated as z/, and K{q, {p,r)) = -^{Vp^{p,r))'^{a{q,r))'^, 
u{p,r) = {v{p,r)Y' in [Zl Lemma 5.1]. We apply this lemma and then we get that for 
all V ^ A° 


'q,u^ 


\v{p,r)\'^p{dpdr\q) > { 1 ^ - rm{q))^Qm{q){i' - r^iq)) 


Jy 

Moreover, if we take 


1 


(23) v{p, r) = u^q, {p, r)) = {q, r)(Vp^>(p, r))'^Qj-{q){iy - r^iq)) , 


m 


we will have 


I 


Wu{q, {P,r))\^ p{dpdr\q) = {v - rm{q)fQm{q){i' - rm{q)) ■ 


This shows that 


= rm{q)VQm\q){i' - rm{q)) , 


and the minimum is achieved in 

Now, that we have identified that the action functional can be written in the pro¬ 
ceeding form we can proceed in proving the Laplace principle upper bound. We must 
show that for all bounded, continuous functions h mapping C([ 0 , T];M'^) into M 

h{q^) 


lim sup —elnE 
e4.0 


exp 


< inf [Sm{(j)) + h{(j))] . 
<^eC([0,T];R‘i) 


By the variational representation formula, it is enough to prove that 


(24) 


lim sup inf E 

£4-0 


1 


jusf ds + h{q^) 


< inf [Smi<P) + h{(p)] . 

0eC([O,T];Kd) 
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To be precise, we consider for the limiting variational problem in the Laplace prin¬ 
ciple a nearly optimal control pair (?/),P). In particular, let t/ > 0 be given and consider 
V’ G C([0, with i/jq = qo such that 


<S'm(V’) + < inf [Sm{(p) + h{(p)]+r] <oo. 

0eC([o,r];K‘^) 

It is clear now that L°{x, v) is continuous and finite at each pair (x, v) G Hence, 
a standard mollification argument, allows us to assume that is piecewise constant, see 
Lemmas 6.5.3 and 6.5.5 in Subsection 6.5 of [6]. The control in feedback from used to 
prove is then given by ([23]) . i.e. 


It is easy to see that Condition 12.11 guarantees that ut is continuous in all of its 
arguments and that d?!) has a unique strong solution with ut = Ut- Then, by Theorem 
13.11 we obtain that in distribution q'^ —)■ q, where 


qt = qo + 



0 Jy 


7 (qs,iP,r),u^^{qs,ip,r))^ n{dpdr\qs)ds. 


Keeping in mind the definition of A° ■ and that ipo = qo, we obtain that 
qt = do + / '>Psds = for any t G [0, T], with probability 1 . 

Jo 

Therefore, we finally obtain that 


lim sup 

etO 


—e InE 


exp 


= lim sup inf E 
etO “ 




< lim sup E 
etO 


1 


2 Jo 


\ut\‘^ dt + h{(f 


= E [5^(X) + h{X)\ 

< inf [Sm{4>) + h{4>)] + r]. 

<i.eC([0,T];R‘^) 

Since q is arbitrary, we are done with the proof of the Laplace principle upper 
bound. At the same time we get the explicit form of the rate function 

(25) 

SM = { 2 I 
+ 00 


- rm{4>s)) Qm {4>s){4>s - rm{(t>s))ds if ()) G AC([0,r];M'^),(/)o = qo 

otherwise . 


Part 3. [Compactness of level sets]. This follows directly from the alternative 
representation (|2^ . as it is in the standard quadratic form, see for example [7]. 


This concludes the proof of the theorem as well as of the alternative representation 
of Theorem 12.11 □ 
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5. Convergence of the action functional as m —0 


Let /3 > 0 and set a{q, r) = ■s/2l3\{q)I. Recall the definitions of the operators 
and and of the corresponding invariant measures from Section [2j 

Theorem 12.21 follows directly from Lemma [5.II below, whose proof is based on The¬ 
orems 12.31 and Eai 

Lemma 5.1. Assume that Conditions [ET\ and lK^ hold. Let Qm{q),rm{Q) cmd Qo{q),rQ{q) 
be as in Theorems \2. 1\ and \2.2\ resvectivelv. Then, for any e > 0, there exist some mo > 0 
such that for every g G and every 0 < m < mo we have 

\Qmiq) - Qo{q)\ < £ , \rm{q) - ?'o(9)| < e • 

Proof. For notational convenience and without loss of generality, we shall set a{q, r) = 
21, fi = X{q) = 1 . Since q is viewed as a parameter, we do not mention it explicitly in 
the formulas. We have 

2 


Qm Qo — 2 


>y 


1 


m 


dp^{p,r) 


p^{p,r)p'^{p,r)dpdr — 2 f \I + drx{f’)f P^iPT'£)dpd: 

Jy 


= 2 


ly 


m 


dp^{p,r] 


- \I + drX{r)Y 


p^{p, r)p^{p, r)dpdr 


+ 2 f \I Adrxir)]"^ {p^{p,r) - 1) p'^{p,r)dpdr 

Jy 


= 2 


ty 


m 


dp<^{p,r) - il + drx{r)) 


m 


dp<^{p,r) + (I + drxir)) 


p^{p, r)dpdr 


+ 2 f \I + drxir)\‘^ {p'^iP,r) - 1) p°ip,r)dpdr . 

Jy 


ty 

Taking absolute value and using Cauchy-Schwarz inequality we obtain 


\Qm — Qo\ <2 


1 


dp^ -{1 + drx) 


1 


+ 2 


L2(y;p0) 


\\P^-nL^i.y;pO) 


dp^ + (/ + drx) 


L^iy-,P^) 


{I + drX? 

By Theorem 12.31 we have 

(26) IIP™-Illy 

By Theorem 12.41 we have 


(27) 


lim 

m—^0 


1 


m 


dp<^ - (/ + drx) 


= 0 


L^(y\P^) 


The results (1261) and (I27p imply that there exists a uniform constant C such that 


(28) 


sup 

me( 0 ,l) 


dp^ 


+ sup 

L' 2 {y-,pm) me(o,i) 


\I + drX\^ 


L^{y;p^) 


<C, 
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and by classical elliptic regularity theory there exists a uniform constant C, clearly 
independent of m, such that 


(29) 


{I + drxf 


< C . 

L^{y-,P°) 


From (I26I) - (I29I) we infer the first inequality of this Lemma. In a similar way from (|26p 
and (|28p we also derive the second estimate of this Lemma. □ 


6. Convergence of the invariant density 

In this section we prove Theorem 12.31 For notational convenience and without loss 
of generality, let us assume in this Section that a{q, r) = 21 and that /3 = A = I (recall 
a{q,r) = ■s/2X{q)pi). Since q is viewed as a parameter, it will not be mentioned 
explicitly. 

We want to show that 


(30) 


hmjr-lhw) 


= 0 . 


where we recall that p^{p,r) = p^^{p)po{r). 

Notice that in the case of gradient potential, i.e., when b{q,r) = —VrV{q,r) then 
(j30p is immediately true even without the limit. In fact in this case we have that the 
invariant density is basically p™'(p, r) = p^^{p)po{r) for every hnite m G ]R_|_ which 
implies that p^{p,r) = 1 completing the proof of (1301) . Our goal here is to show that 
this true in the more general setting of not potential drifts. 

By Condition 12.II the drift b{q,r) and its partial derivatives are uniformly bounded 
with respect to q. For this reason we sometimes suppress the dependence on q and write 
b{q,r) = b{r). Also, for notational convenience, let us set 


h{r) = b{r) - Vrlogpo{r) . 

This definition for h{r) will also be used throughout the rest of the paper. 

Notice that in the gradient case, i.e, when b{r) = — Vy(r), we have that h{r) = 0, 
but in the general case one has h{r) ^ 0. Let us next establish some useful relations 


Lemma 6.1. Let f,g be two functions that belong in the domain of definition of C 
Then, we have the identity 


'y 


[(^™/(p, ?’)) gip, r) + {C^gip, r)) f{p, r)] p^{p, r)dpdr = 


- — f ^pf{p,r)Vpg{p,r)p^{p,r)dpdr + -^ [ f{p,r)g{p,r)h{r)pp°{p,r)dpdT 
m Jy ^ y/m Jy 

In particular, we have that 

[ r)) f{p, r)p^{p, r)dpdr = 

Jy 

= -— [ \Vpf{p,r)\'^p°{p,r)dpdr + -^ [ \f{p,r)\'^h{r)pp^{p,r)dpdr. 
fn Jy Zy/m Jy 


19 









Lemma 6.2. Let f,g be two functions that are in yv\’^{y), i.e., the functions and their 
derivatives with respect to p are in L'^{y). Then, there exists a finite constant K < oo 
that depends only on sup^g^'^ l^(^)l such that 


{h{r)p, fg)L 2 (y.pO 


< K 


ll/llL2(:y;pO) \\'^p9\\L^y-,pO) + W'^pfWi^yipO) \\9\\L^(y-pO 


Lemma 6.3. For every g > 0, there exists constant constant K < oo that depends only 
on l^(^)l such that 


1 


> -K 


1 


9 \\f\\Lfiy-,pO) + ^ W'^p/hfiy-po 


(/>^/)L2(y;pO) = 2 \f\^)L^{y-,pO) 

where we recall that Bf = p ■ Vrf + b{q, r)Vpf. 

The proof of Lemmas I6.1H6.3I are in Appendix B. Let us now define 

<5-(p,r) = p-(p,r)-l. 

Recall that our goal is to prove Theorem l2.3l i.e. that (I30p holds. The next lemmas 
are towards this direction. The proof of Lemmas l6.4H6.6l are in Appendix B. 

Lemma 6.4. For every m > 0 we have the following equality 

,2 ^An 


IV 

ApO ||^2(y.p0) 


(h(r)p, + yM {h{r)p, S^LW) • 


L^y-,pO) 

Lemma 6.5. There is a universal constant K > 0 that depends on sup^g-j-d \h{r)\, but 
not on m > 0, such that for all m sufficiently small 

(1 - V^) ||VpVp5-||i2(y;,0) < V^K [l + + \\^pSn\l^y.,pO) • 

Lemma 6 .6. There is a universal constant K > 0 that depends on sup^g-j-d max(|/i(r)|, |Vr/i(r)|), 
but not on m > 0, such that for all m sufficiently small 

(1 — \/m) ||VpVr( 5 ™'||p^ 2 (y.pO) ^ 

< ^/mK 1 + ||(^”'||i2(;y.p0) + l|Vp5”'||^2(;y.pO) + 11 5"* 11 ^2 (y .pO) • 

Let us define to be the operator with m = 1. We recall that 

C^=A + B, 

where A = —p ■ Vp + Ap and B = p ■ Vr + b{q, r) ■ Vp. It is easy to check that, with 
respect to the measure p^{p,r)dpdr we can actually write that 

= -AA* + B 

where 

A = Vp , and A* = — (Vp — p) . 

One can also check that the adjoint operator of B is formally given by 

* = —B + ph(r) . 
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Notice that the latter relation implies that B is antisymmetric only if h{r) = 0 
which essentially is the case of gradient drift. However, in the general case h(r) ^ 0 
which would imply that B is not antisymmetric. Next, we introduce the operator 


C = [A,B] = [Vp,pVr + 6(r)Vp] = . 

A word on notation now. In order to make the notation lighter we will write from 
now on 

INI = IMIl2(:V;p°) ’ ^ ■)L2(y;p°) ’ 

for the norm and for the inner product in the space L?‘{y; p^). 

In order to show that ( 1301 ) holds, we use the work of [ 30 ]. In particular, as in m, 
let a, 6 , c be constants to be chosen such that l>a>b>c>0 and let us define the 
norm 


((/, /)) = ll/f + « \\Aff + 267^ {Af,Cf) + c \\Cff . 


In fact, as it is argued in m, the norms ((/,/)) and ||/||^i(y.pO) are equivalent as 
soon as 6 < y/ac, in that 

min{I, a, c} (^1 - Il/llb(y;p 0 ) < ((/, /)) < max{I, a, c} (^1 + Il/llb(y;p 0 ) • 

Since, we are dealing with a real Hilbert space, all the inner products are real. By 
polarization we have 


((/,£V)) = {f,C^f) + a{Af,AC^f) + b [{AC^f,Cf) + {Af,CC^f)] +c{Cf,CC^f) 


One important difference between the current setup and the setup of |30j is that 
there B* = —B whereas here that is not the case, as we have B* = —B + ph{r). Keeping 
that in mind and repeating the argument of the proof of Theorem 18 in [30], we obtain 
that there are constants a, b, c that are sufficiently small such that 1 a S> 1 ) S> 2 c with 
6 ^ < ac (the exact same constants as in [30]) such that 

{{f,C^f))>K[\\Aff + \\Cff 

+ {(/, Bf) + a (H/, BA!) + b {phir)Af, Cf) + c {Cf, BCf)} 


(31) 


> ^ II/IIhi + {(/, Bf) + a {Af, BAf) + b {ph{r), AfCf) + c {Cf, BCf)} . 


The bracket term of the right hand side of the inequality is due to the fact that in 
our case h{r) 0 and thus B is not anti-symmetric. The bracket term is equal to zero 
in [30]. 

Let us now choose / = d™' in (I31|) . The strategy of the proof is: (a) bound from 
below the bracket term on the right hand side of (1311) using Lemmas 16.2116.61 and the 
equation that satisfies, and (b) bound from above the left hand side of (l3T]) using 
Lemmas 16.2116.61 and the equation that <5™ satisfies. Putting the two bounds together 
one will then obtain a bound for which will give the convergence to zero of (130^ 

that we need, combined with Poincare inequality for the measure p^{p,r)dpdr. 
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We would like to highlight here that one of the obstacles in putting the lower and 
upper bounds together, are the order one terms (/, in the definition of 
and {f,Bf) in the lower bound (l3T]) . However, as it turns out, see (1351) . for / = 5™, we 
actually have that (£^5™, 5™) — {BS'^,5'^) = o{y/m) which then allows us to proceed 
with the bounds. The rest of the terms are being handled via Lemmas I6.2ti6.6l 

We start with obtaining a lower bound for the bracket term on the right hand side 
of ([311) using Lemmas I6.4H6.6I and the equation that 6'^ satisfies. For this purpose, let 
us define 


^ +a{A6'^,BA6"^)+ b{ph{r)A6^,a"')+c{C6'^,BC6^) 

= {6^,B5^) + Ri{S^) . 

Let 77 > 0 to be chosen. By Lemmas 16.2116.31 recalling that Ad"^ = VpJ™ and 
C5'^ = and using the generalized Cauchy inequality ab < pa? + ^ 6 ^ we have that 


>-K{a 
+b 


’)IIW'V‘ll" + L||v,v,^” 


+ b 


\Vr5^\? + ^||VpVp5^ 


4?7 


1 


ri\\Vp6^f + —\\VpVr5'^f 


+ c 


1 


V\\Vrd^f + —\\VpVr5'^f 


v\\Vp6"^f + ^ (1 + + ||Vp,5”^f) 


p\\Vr5^f + ^ (1 + + l|Vp<5"*|n 

4?7 / 


p\\Vp6^f + ^ + l|Vp<5™f + ) 


Next, using Lemmas 16.4116.51 and 16.61 we subsequently obtain 
Ri{S^) > |a 
+b 
+b 
+c 

> -K 

> -K 


77||V,(5™f + ^ + ||Vp5”^f + ||V,<5"^|n 

Arj ^ 


1 (l|v/’"||^ + ||Vr«’"||^) + ^ (1 + + l|V,i”|p + l|VrS”f) 


nn\h + ??{i + \\nf + \\n\h) 

where the positive constant K < 00 may change from line to line but it is always 
independent of m. Choosing now p = p{m) such that hmm 4 .o p{m) = lirnm^o = 0, 
we obtain for 77 ( 771 ) = max{ 77 ( 77 i), 7 ^^} i 0 , that 

Ri{6n > -Kp{m) [1 + + lld-ll^i] . 

So, overall we have that for m sufficiently small there is 77 ( 771 ) 0 as m 4- 0 such 

that 

R{5^) > {5'^,B5^) - Kp{m) [l + + \\5^\\jpr] , 
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or in other words by (I3ip with / = 5”^ we have that for m sufficiently small there is 
fj^m) I 0 as m I 0 such that 

(32) {i6^,C^6n)>K{l-fi{m))\\5^\\l^ + {5^,B6^)-Ki)im)[l + \\6"^\\^] . 

Hence, recalling the dehnition of the inner product ((•, •)), using (I32p and rearrang¬ 
ing the expression a little bit we have obtained the following bound 

K(1 -r}(m))||5™||^i < (5™,£^<5”^) - [l + ||(5™||2] 

(33) + + +c(C<5”^,C£M”^) . 

The next goal is to derive an appropriate upper bound for the left hand side 
of (p3l) . First, we need to obtain the equation that satisfies. By factoring out 
p^{p,r) = p°(p, r)p™'(p, r) where p^{p,r) = p^^{p)po{r), we obtain the following equa¬ 
tion for p^{p,r): 

£^nP,r) = -^Bp^{p,r) - -^ph{r)p^{p,r) . 

^ y/m y/m 

where we recall that h{r) = b{r) — Vrlogpo(^)- Hence, the equation for d^{p,r) = 
p^{p, r) — 1 is 

£™5"^(p,r) = -^BS^{p,r) - -^ph{r) [6^{p,r) + 1] , 

^ y/m y/m 

or in terms of the operator we have 

(34) C^6^{p,r) = (l + y/m) BS^{p,r) — y/mph{r) [6"^(p,r) + 1] . 

By multiplying both sides of (1341) by <5™ and integrating over y with respect to the 
measure p^{p,r)dpdr we then obtain that 

(35) {C^5^, 6^) - {B6^, 6^) = y/^ {B6^, 6^) - {ph{r), {S^ + 1) 6^) . 

Hence, using (1331) and (|35P we have the following bound 

A'(l - r)(m)) ||<5”^||^i < {B6^, 6^) - {ph{r), {6^ + 1) 5^) + Kfi{m) [l + ] 

+ a{A5'^,AC^5^) + h[{AC^5^,C5^) + {A5^,CC^5^)] A c{C5^,CC^5^) 

<Kf,{m) [l + ||<5”"f] 

+ [{B5^, 6^) - (ph(r), (5™ + 1) 5^)] 

+ a{A5^,AC^5^) + b[{AC^5'^,C5^) + {A5^,CC^5^)] A c{C5'^,CC^5^) 

<Kf,{m) [l + ||5”"f] 

(36) 

+ Ti{5^) + aT2{6^) + bTsib^) + 0^(6^) . 

Our next goal is to derive upper bounds for the terms for z = 1, 2,3,4. For 

better readability, we collect the required bounds in the following lemma, which we also 
prove in Appendix B. 
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Lemma 6.7. Let the terms Ti{5'^) for i = 1,2,3,4 be defined as in i36\) . Then, there 
exists a constant K < oo that does not depend on m, and a sequence ri{m), 0 as 

m —0 such that for m sufficiently small the following bounds holds 


\T,{5n\<v^^\\n\ 

| 72 ( 5 ™)| < ^{m) + ^/m + 
| 73 (<^”*)| < (^ri{m) + y/m + 



Now that we have obtained the desired bounds for the terms for i = 1, 2, 3,4 

let us put them together. There are some constants Ki,K 2 < oo, and a sequence 
f}{m) = max{r7(7T7-), i 0 such that for m sufficiently small 


(1 - 7}(m)) < fi{m)Ki [l + ||(5™||2] + \^Ki + y/ffiKi ||Vp<5" 


+ a 
+ b 
+ c 


rm||2 


fj{m)Ki\\6"^\\jji + fi{m)Ki (1 + \\5' 


fiim)Ki i + p-f + +{l + y/fil)\\Vr6^ 




(37) 


< (^(m) + V^)i72 l + ||5’"|r +115"* 11^1 + 6(1 + ||V,5" 


Now we choose m small enough such that fj{m) < 1, {fi{m) + y/rn)K 2 < 1/2. 
Moreover, we also note that since by construction 6 <C 1 we can write for m small 
enough 6(1 + y/m) <C 1/2. In fact the proof of [30] shows that we can choose a, 6, c to be 
positive but as small as we want, as long we choose the constants a, 6, c to be ordered 
appropriately. Putting these estimates together, we get that there is some constant 
it's < oo such that for m small enough, one has 


(38) 


l|5^ 


[i + r 


1/2 — fi{m) 

In order now to close the estimate we need to use Poincare inequality. Here we 
make the assumption that the drift 6(r) is such that the invariant measure p^{p,r)dpdr 
satisfies the Poincare inequality with constant k > 0 . In particular, for a function 
Q{p,r), we have that the Poincare inequality in the following form holds 

2 


Q- / Q 

Jy 


L^{y-,P°) 




Let us set now Q{p,r) = 6"^{p,r). Notice that by definition of 6^{p,r) we have 



d^{p,r)p^{p,r)dpdr = 0 . 
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Therefore, we have obtained 


(39) 




12 ^ 11 rm 112 

I 11^1 


Inserting now (1391) into (1381) . we finally obtain that for m small enough 

(40) , 

from which the desired result finally follows: 

pm ||2 ^ {fl{m) + y/m) —)■ 0 . 

This concludes the L‘^{y-,p^) convergence of the invariant measures. 

7. Convergence of the solution to the cell problem 

The goal of this section is to analyze the cell problem ([8]) that <I>(p, r) satisfies and 
we want to prove Theorem 12.41 As it will become clear from the proof below, we prove 
even more. We rigorously derive an asymptotic expansion of $(p, r) in terms of powers 
of y/m. 

Let us recall our assumption a{q,r) = 2j3\{q)I. Let I = l,2,...,(i be a given 
direction and let us define 


the equation that 


Tfl 

r) = d>£(p, r) - ■ et , 

where ei is the unit vector in direction 1. Then, bearing in mind 
y£{p,r) satisfies is given by 

(41) £”-P,(p.r) = -l^ . 

Moreover, by Condition 12.21 we have that that for every m > 0 

y£{p,r)fj,{drdp\q) = 0 . 

Let us write for notational convenience the hypoelliptic operator 


L 


£T = ^A+^B , 


^ m yjm 

where we have already defined A = —p ■ Vp + /3Ap and B = p - Vr + b{q, r) ■ Vp . 

Let us now write the expansion 

r) = r) + y/my£p{p, r) + r) + d'™ 3 (p, r) . 

Assume that the functions ^and satisfy the following equations 

(42) Ayefi{p,r) = 0, 

(43) B^ifiip, r) + X{q)Aye,i{p, r) = 0 , 

(44) B^ipip, r) + X{q)Aye,2{P^ d = -biid, r) , 

(45) C^-^ddP^ r) = - V^^T,, 2 (p, r) . 
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and that for i = 0,1,2,3, we have that / ^ii{p,r)p^{p,r\q) = 0. The next step is 

Jy 

to analyze the solutions to ()l2]) - (H5]) . First we notice that (H2]l basically implies that 
i-e., function is a function of r alone. Then, using this we 

get by (fl3l) that it has to be the case that 

r) = o(r-) • p + ^£,o{r) 

Kq) 


for some function 'l/^o(^)- From equation (HIT) and (H2]) . (H5t) . (03]), (H5|) we see that 
up to an additive constant we can assume that T£ o(^) = 0. Lastly, we notice that the 
solvability condition for (|44l) is 


/ 

JR” 


[B^e^i{p, r) + be{q, r)] 7r{dp) = 0 ^ 

B{Vr'^£fi{r) ■ p) + B'^£fl{r) + b£{q, r) 7r{dp) = 0 ^ 

+ HQ^r) ■ Vr'^£^o{r) + p-Vr'^£fi{r) + b£{q,r) TT{dp) = 0 


l3Ar'i>£fl{r) + b{q, r) ■ Vr^£fi{r) = -bi{q, r) , 


.kr 


where the Gaussian structure of the invariant measure Ti{dp) ~ e Ts” dp and integration 
by parts were used. Notice that this is exactly the solution to ([9]) with a = 2j3X{q)I. 
Thus, by uniqueness of the solution to dOj) we basically have that for every i = 1, - ■ ■ ,d 


^£,o{^) = X£ir)- 


Hence, we have that 


Vp4>(p,r) = Vp 

= V. ( 


Kq) 


p + 'I'(p, r) 


p 

m 


p + To(p, r) + ^/m^i{p,r) + mT 2 (p, r) + r) 


K<l) 


KKaY 

[I + VrXir)] + mVp4'2(p, r) + VpT™(p, r) 


Having established the last display, it is easy to see that in order to show (l27)l . we 
basically need to show that 


lim 

m —>0 


VmVp'i/2{p,r) + ^VpT^(p,r) 

-v/m 


= 0 


L^{y-,p^) 


or, in other words, it is sufficient to show 

(46) hmJ|VTOVp4'2(p,r-)||p^2(y.^^) = 0 , 

and 

1 


(47) 


lim 

m—^O 


Vp4/^(p,r: 


= 0 


L2(k;p-) 
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Relation ([l 6 t) can be claimed to be true by the fact that r) is solution to the 

elliptic problem (I44p and Theorem 12.31 

So, it remains to prove (|T7l) . At this point let us recall that ^ 73 (p, r”) is solution to 
(|15|) . i.e., it solves 


(48) 




Notice that for the purposes of this section q is seen as a fixed parameter by the 
operators and recall that we have already assumed a{q, r) = 2j3\{q)I. Namely /3X{q) is 
seen as a fixed constant. Hence, from now on and for notational convenience, we shall 
assume without loss of generality that a(q,r) = 21, i.e., that P = X{q) = 1. Let us first 
apply Lemma 16.11 and we get 


(49) 

Jy 




1 


m 




A 


{'^e^3?Hr)pp^{p,r)dpdr 


Lemmas 17.1117.41 that follow are proven in Appendix ICl 


Lemma 7.1. We have the Poincare inequality 


(50) 


/- 


/ f{P,r)p^{p,r)dpdr 

Jy 


< K IIV 


L^(y-,pO) 


pf\\L^{y-pO) 


for some constant k> 0 independent ofm. 


Lemma 7.2. We have 


(51 j hmsup-- < O < oo ,hmsup 

m—>0 ITT' ' m—>0 


V 4/™ 

^p^e,3 


L\y-,p°) ^ ^ 


< oo 


m 


for some constant C > Q independent ofm. 


Lemma 7.3. We have 


(52) hmJ|4'^3||i4(y.^o) =0 . 

Lemma 7.4. We have 

{J^'^f{p,r))g{p,r)p'^{p,r)dpdr+ I f{p,r) g{jp,r)) p^{p,r)dpdr = 

jy 

= I [Vp/(p, r) • a{q, r)Vpg{p, r)] p"^{p, r)dpdr 

mjy 

We set in particularly in ()49p / = g = '£'^ 3 , then we will have that 

[ ^e^3dr)p^{p,r)dpdr =-— [ \Vp^^^{p,rdp^{p,r)dpdr . 

Jy m Jy 


Jy 

(53) 
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But, we also know that 'if'^^{p,r) satisfies (HHI) . Therefore, multiplying both sides 
of (j48p by r) and integrating against the invariant density p'^{p,r) gives us the 

identity 

— [ |Vp^'™3(p,r)|2p”"(p,r)dp(ir = / {B'^i^2{p,r))'^'l^3{p,r)p'^{p,r)dpdr , 

m Jy 

or, in other words 
(54) 


>y 




L2(y;p"*) 






We now have the estimate 


I ^£,3 


Ih^y.pm) = j^{^^3fp"^{p,r)dpdr 

= [ i'i' 7 : 3 ?^"^iP,r)p\p,r)dpdr+ f p\p,r)dpdr 

Jy Jy 

< ll'J'b3lli^(y;pO)l|5”^lli.(y;pO) + • 


Applying Lemma [TIBI and Lemma \T?I\ and the fact that lim \\d'^\\L'^{y-pO\ = 0 we 

m—>-0 ' 


see that 




Thus we have by 


1 


Vp4'^3 


< 


L2(y;p’") 




as m ^ 0. This is ()47p . completing the proof of Theorem 12.41 


Appendix A. On properties of the solution to the hypoelliptic cell 

PROBLEM 

In this section we recall some results on the solution to the hypoelliptic Poisson 
equation ([SP from [T3]. Since the set-up of the current paper has some differences from 
the setup in m, we formulate the results that we need in the current setup, even though 
we emphasize that the derivation follows basically from mi- 

under the assumptions made in this paper, Theorem 3.3 from m guarantees that, 
(jSP has a smooth solution that does not grow too fast at infinity. In particular, we have 
that for every p > 0, we can write 

<I>(p, r) = e2lPl^$(p, r) 

where € S, the Schwartz space of smooth functions with fast decay. Further¬ 
more, as it can be derived from the proof of Theorem 3.3 of mi, if we let Umax = 
maxjj=i^...rf sup(qr)|, then we have that for every p G (0, 2(T“ax) f^ie solution <I> 
is unique (up to additive constants) in the space (y, dpdr^. 
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Moreover, it is clear that for each fixed g, the operator defines a hypoelliptic 
diffusion process on (p,r) S 3^ = x T'^. Let us define this process by {p.,r.). We 
recall then the following useful bounds from |14j . 

(i) There exists a constant C such that 


E 


e 2 



< E 


-2 

^max 1^ |2 


Ipof+Ct 


(ii) For every T > 0, there exist constants r],C > 0 such that 


E 

sup 

< CE 



te[o,r] 




Based then on these bounds, the computations of m reveal that the following 
bounds for the solution to (l8|) are true. In particular we have that for every T,p > 0 there 




;£ 

<5 


< 


exists a constant C > 0 that is independent of e, 6 such that E sup^gjQ 

and E sup^gjo^T] Vp4> (pf, ^ <C. 

These bounds are used in the proofs of this paper and more specifically in the 
derivation of Theorems [Q and 14.11 


Appendix B. Proofs of Lemmas in Section [3 

Proof of Lemma \6.1[ The proof goes in a standard way using integration by parts. We 
present the main steps for completeness. 


ly 


?’)) 9{P, r)p^{p, r)dpdr = / f{p, r) [g{p, r)p^{p, r)) dpdr 


'y 


dpdr 


= f f{p,r) ^{A)* {p°{p,r)) g{p,r) + {-C^g{p,r)) p^{p,r) 

+ ^ j^f{p,r) [(p°(p,r)/ : Vlg{p,r) + Vpp°(p,r) • IVpg{p,r))] dpdr 

-E f r)g{p, r) [6(r)Vpp°(p, r) + pVrp°(p, r)] dpdr 

y/m Jy^ 


ly 

= / f{p,r){-CTg{p,r))p°{p,r)dpdr 

Jy 

H-/ fiP, r) (P°(P; r)I ■ Vp 5 (p, r) + Vpp°{p, r) ■ IVpg{p, r)) dpdr 


m 


+ 


jy 

-j= f{p,r)gip,r)p-h{r)p^{p,r)dpdr 

Am Jy 


= / f{p,r){-Cf^g{p,r)) p'^{p,r)dpdr 


ly 


- — I P) • I'^pdiP, r)] P^iP, r)dpdr + / f{p, r)g{p, r)p ■ h{r)p^{p, r)dpdr . 

m Jy Jy 

To derive the last line, we used integration by parts as well as the definition h{r) = 
b{r) — Vr log poA)- The statement of the lemma follows. □ 
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Proof of Lemma \6.SX We start with the following calculation 


V, 


(e =-pe f{p,r)g{p,r) + e (/(p, r) 5 r(p, r)) 


Therefore, we obtain 

/ Pf{p,r)g{p,r)p^^{p)dp= [ Vp{f{p,r)g{p,r))p^^{p)dp 

= / (Vp/(p, r)g{p, r) + f{p, r)Vpg{p, r)) p^^ip)dp . 

jRd, 

Multiplying both sides by h(r)pQ[r) and integrating over r S we then obtain 
after using Holder inequality 

{Kr)p, fg)L^{y-pO) = (Vp/(p, r)g{p, r) + f{p, r)Vpg{p, r)) p°(p, r)dpdr 


< K 


Il/lli2(y;p0) l|Vp5llL2(y.pO) + l|Vp/||^2(y.pO) 11511^2(^.^0) 


This completes the statement of the lemma. 
Proof of Lemma [gTM We notice that 


□ 


{f,Bf)L^{y-pO) = J^pVrf{p,r)f{p,r)p^{p,r)dpdr 

+ / b{r)Vpf{p,r)f{p,r)p°{p,r)dpdr 

Jy 

= Termlrn + Term2m ■ 


By integration by parts, we have 
Terrnlm = / pVrf{p,r)f{p,r)p^{p,r)dpdr 

Jy 


= - pf{p,r)Vrf{p,r)p^{p,r)dpdr- / pV^ logpo(?’) |/(pt)I p ip,r)dpdr 

Jy Jy 


Thus, we get 


1 


Terrnlm = J^pVrf{p,r)f{p,r)p'^{p,r)dpdr = J^pVrlog po{r) \f{p,r)f p^{p,r)dpdr . 


Similarly, we have 

Term2m = [ b{r)Vpf{p,r)f{p,r)p^{p,r)dpdr 

Jy 


= — J b(r)/(p,r)Vp/(p,r)p {p,r)dpdr + j fc(r)p|/(p,r)| p {p,r)dpdr 


>y 


Thus, we get 


1 


Term2m = / b{r)Vpf{p,r)f{p,r)p^{p,r)dpdr = - / 6(r)p|/(p,r)| p'^(p,r)(ipfir . 

Jy ^ Jy 
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Putting the representations of Termlm and Term2m together, we have in fact ob¬ 
tained 

(/,^/)L2(y.pO) = ^ J^pHr)\f{p,r)f p^{p,r)dpdr . 

Hence, by Lemma 16.21 we have that there exists a constant K < oo that depends 
on sup^g^d \h{r)\ such that 


(/)^/)L2(y;pO) — 2 l/P)L2(y;pO) > ll/llL2(y;pO) 11 ^p/ 11 ^2 (;y .pO 


> -K 


d\\f\\LHy-,pO) + ^ l|Vp/llL2(y;pO) 


where we use the generalized Cauchy-Schwarz inequality ab < for any 

r/ G (0, oo). This concludes the proof of the lemma. □ 

Proof of Lemma \6.4\ Recall that by ([6|), the equation for d'^(p,r) = p^{p,r) — 1 is 

(55) r) = -^B5^{p, r) - -^ph{r) [6^{p, r) + 1] . 

^ \/m wm 


Let us multiply now the last equation by S"^{p,r) and integrate over y against 
p^{p,r). Doing so, we get 

^ (Bs^.nLW) - (?'"(’■) ■ 

The next step is to rewrite the term Lemma [6l^ we have 

(^/>/)L2(y.pO) = ^ J^ph(r)\f(p,r)\^ p°(p,r)dpdr . 

Inserting the latter expression into (1561) we obtain 


1 


{ph(r) lf(p,r) + 1 ] ,/)^2(y.p0) 


<PMO,l/P>L2(y.p0) 


(ph(r) [f(p,r) + 1 ] ,/)L2(;y.pO) 


{pKr)J)L^(y-,pO) ■ 

Next step is to apply Lemma [6T] with f{p,r) = g{p,r) = 5 ^{p,r) to get 

Combining the last two expressions, we obtain 


(56) - 


1 


m 


I 11^ I 

llL2(:y;pO) + 


2^/m 


h{T)p,\5^\^) 


L^iy-,pO) 


{Ph{r),6^)^2^y.p0) 
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and after rearranging, we obtain 


I Yy cm 112 

I \\ L ^{ y ; pO ) 


h{r)p, 0 ) + ^ {Ph{r), ^'^)L^y-,pO) 


2 \ " ' /L 2 {y-,p 0 ) 

This concludes the proof of the lemma. 


□ 


Proof of Lemma |g.51 The proof goes along the same lines of Lemma 16.41 We take ^p^ 
on both sides of the equation ([55]) and we get the following equation 


1 


1 


CLdp^5^ = - —hiir)[5^ + 1 ] - • h{r)dp^S 


/m y/m y/m 

Multiplying both sides of the above equation by dp^ 6 "^ and integrate with respect 
to L‘^{y;p^ )-inner product we get 

{C^dp^5-^,dp^5-^)L2^y.^pO) = 


m 

1 


{hi{r){S^ + l),dp^ 6 nLHy-,Po) 

{p-h{r)dp, 6 ^,dp^ 6 ^)L 2 ^y.p^) . 


We apply Lemma lUTTl with f{p,r) = g{p,r) = ^p^5'^{p,r) to get 

{£^a,,s’-,d,,nLW> = . 


We now apply Lemma 16.,41 and we have 

{Bdp^6'^,dp^6"^)L2(^y.p0) = ^ J^p- h{r)\dp^ 5 ^\^p°{p,r)dpdr . 

Furthermore, we can calculate 

{h,{r){d^ + l),dp^6nL2^y,p0) = [ h,{r){6^ + l)dp^S^p\p,r)dpdr , 

jy 

{p ■ h{r)dp^6'^, dp^5'^)L2(^y.pO) = [ p- /i(r)| 0 p.< 5 ""|V°(p, r)dpdr . 

Jy 

Thus, we get the identity 


1 


2 y/m 


{ph{r), \dp. 6 ^\^)L 2 (y-p 0 ) 
{h,{r)i6^ + l),dp^5^)L2^y,pO) . 


__||V /) 

^W'^pOpiO \\L2{y-pO) 


Making use of Lemma 16.21 and Young’s inequality we estimate 

+I<Ai(\\A\Uy.,y.) + + 1) . 

where iL > 0 is a constant that depends only on sup \h{r)\. This implies the lemma. □ 
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Proof of Lemma \6.(A The proof goes again along the same lines of Lemma 16.41 We take 
dr^ on both sides of the equation ([55]) and we get the following equation 




^ Bdr.8'^ - ip • + 1] - ip • h{r)dr^5^ 


m 


m 


m 


Multiplying both sides of the above equation by dr^d^ and integrate with respect 
to L^{y;p^ )-inner product we get 

^ -.{p-drMr)[S^ + M,dr^Snmy-,po) 


m 


1 




We apply Lemma IH tI with /(p, r) = g{p,r) = dr^5'^{p,r) to get 

1 

2y/rn 


{ C ^ dr , d ^, dr ,5^) L 2^ y ., pO ) = . 


We now apply Lemma 16.31 and we also have 

{Bdri6^,dri5^)L2(^y.pO) = ^ j^p ■ h{r)\driS'^\‘^p°{p, r)dpdr . 

Furthermore, we can calculate 


(p • (9r,h(r)[5™ + l],driS'^)L 2 (y-pO) = f p- drih{r){5'^ + l)dr ^5^^ p^{p,r)dpdr , 

Jy 

We can apply a straightforward generalization of Lemma 16.21 with h{r) replaced 
by drM{r), as well as Young’s inequality, to estimate the right hand side of the above 
equation by 


{p-drMr)[S^ + ^],dnnLW)^ 


7^ / 11 rm 112 _i_ 11V7 112 11V7 A”^ 112 

I II® llL2(y;pO) + II VpO llL2(y;pO) IIllL2(y;pO) 

1 

-L—IIV rl 

+ 211^^*^^*® llL2(y;pO) > 


+ 1 


where iL > 0 is a constant that depends only on sup \Vrh{r)\. We also have 

rST*^ 


(p • h{r)dri5'^,dri5^)L2(y-pO) = j p- h{r)\dri5'^\‘^ p'^{p,r)dpdr . 

Thus we get the identity 

-^llW^n<^™lli2(y;pO) = -^^(ph(r),|9,,(5"*|2)p^2(y.p0) 

-^(p • drMr)[5^ + 1], 9r.<5”^)L2(y;pO) . 

Making use of Lemma 16.21 and Young’s inequality again we estimate 

^pdr,^^\\h^y;pO) < K^\\Vp5^\\l,^y.^p,^ 

+ K^/m{\\5'^\\\2^y.pO^) + IVp'^™llL2(y;pO) + IIW(^™||^2(3;.pO) + 1) 
+ Vm\\V pdr ^ d '^ fp ^^^ y . pO ) , 











where X > 0 is a constant that depends only on sup max(|/i(r)|, |Vr/i(r)|). This implies 


the lemma. 


□ 


Proof of Lemma 6.1. We start with By Lemma 16.31 with f = 6 "^ we have 


T,{5^) = 5^) - {ph{r), {5^ + 1) 5”^)] 


-(ph(r),|5-|2)-(pL(r),(5- + l)n 


= \/m 


= \/m 


--{ph{r), |(5”"p) + {ph{r), 


Thus, by Lemma 16.21 with f = g = 5^ we have the following bound 


|ri(5™)| < V^K ||5”^|| ||Vp(5™|| + ||5^ 


(57) 


-K + 2 


-K 


< y/m —^— WS'^W + Vm— ||Vp(5' 


m ||2 


Next we derive an upper bound for T 2 { 6 ^) = (^Ad"^,A£^ 6 ^'). For this purpose we 
first notice that 


{A6^ , AB5^) = {Vpd^ , VpB5^) 

= {Vp5^,Vp (pVrS^ + b{r)Vp6^)) 

= {Vp6^,BVp6^) + {Vp6^,Vr6^) 

= i {ph{r), lVp5-|2) + {Vp 6 ^,Vr 6 ^) , 

where in the last inequality we used Lemma 16.31 Then, using the equation for 5™, ()35p 
and Lemma 16.21 we have 


|r2(<5”^)| = \{A6^,A£^6^)\ 

= {phir), |Vp(f™|2> + (1 + V^) {Vp5^, VrSn - Hr){6^ + 1)) 

-y/m{ph{r), 


{ph{r), |Vp<5™|2) + (1 + (VpJ™, V,n - ^/^ {Vp 6 ^, h{r){ 6 ^ + 1)) 
< ^^—^^K\\VJ^\\\\VpVp 6 ^\\ + + r/(l + y/f^)\\Vrd^f 


4r/ 


+ {\\d^f + WVpS^f) 




1 


|Vp<5-f+ ^||VpVpJ-|| 1 + 


1 + y/m 
4g 


|Vp<5”^f 


+ r/(l + V^)\\Vr 6 "^f + y/fflK (||,5”^f + ||Vp5™||2) . 
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Next step now is to use Lemma 16.51 Doing so we get the bound 


\T2{5'^)\ = \{A5'^,AC^5^)\ 

(1 - y/m) 


< 


< 


-K 


m 


V/ II + 4 , 


K 


l + ||5”^f+ ||Vp5^ 


■mil 2 


, l + V^IIV7 ^m||2 

+^;r" '■ " 


+ r?(l + V^)\\VrS'^f + V^K{\\5"^f + ||Vp5™f) 

+ i±^'l \\v,s-‘f+ 


8r/ 


Arj J 


+ 7/(1 + V^)||V,<5"*||2 + iL ^ (^1 + . 


Use now Lemma 16.41 and then Lemma 16.21 to bound the term ||Vph’"|p by 

terms of the form y/m (||5™|P + ||Vph™'|p). Choosing then rj = 'q{m) such that 

r]{m) —>• 0 and —)■ 0, we get that for m sufficiently small 

\T2{6^)\ = \{A5^,ACH^)\ 

(58) <!,(„)+ VS;+^)/r||5™||/,. + (^ + V5s)K(l + ll«’"f) . 

Next we derive an upper bound for Ti^{5^) = For this purpose we 

first notice that 

{C6^,CI36^) = {Vr6^,VrB6^) 

= {VrS^yr {pS 7 rS"" + 6 (r)Vp( 5 ™)) 

= {Vr6^,BVrd^) + {Vrd^,Vrb{r)Vp6^) 

= I {ph{r). |V.<5"^|2) + {Vr6^, Vrb{r)VpS^) , 

where in the last inequality we used Lemma 16.31 Then, using the equation for J™" 
T4{6^) = {CS^,CB^S^) 

= {ph{r),\Vrd^\^) + (1 + {Vr5^, Vrb{r)VpSn 

- y/^ {Vr6^,pVrh{r){6^ + 1)) -V^{ph{r), |V,.5”^|2) 

= {ph{r),\Vr6^\^) + (1 + V^) {Vr5^, V,6(r)Vp5™) 

-V^{Vr6"^,pVrh{r){5^ + l)) . 
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Using Lemma 16.21 we subsequently obtain 


\n{5n\ < K\\Vr5^\\\\S/pVr6^\\ + (1 + V^) \{Vr6^,Vrh{r)Vp5^)\ 

+ ^KV,(5™,pV,/i(r)(5™ + 1))| 


< 


(1 - y/m) 


1 


■m\\2 


+ (1 + ^/m) 
(1 - ^/m) 


< 


K\^r^\\Vr5'^r + —\\VpVrS- 

V \\Vp6^f +V^\{VrS^,pVrhir)i5^ + 1))| 

4:7] 

1 


K [v\Nr6"^f + —\\VpVr6^f 

v\\Vr6^f + ^K\\Vp6^f 

+ V^K [||V.5”^|| ||Vp5”^|| + ||VpV.<5”‘|| ||,5'"|| + ' 

(1 - i/m) 


< 


^-K (^rjWVrS^f + j^WVpVrS^f'^ 


r]\\Vr6^f + ^K\\Vp6^f 


+ ^/mK 

(1 - ^/m) 


^ ^ ||Vp<5”^f + ^ ||VpV,5’"f + ^ + ||V,<5^ 


mil 2 


< 


K (^rj\\Vr6^f + ^\\VpVr6^f'^ 


+ (1 + ^/m) 


v\\VrSn\^ + -K\\Vp6' 


+ V^K ||V,<5"^f + ||Vp5™f + ||VpV^<5”^f + ||<5" 


The constant K may change from line to line, but it is always independent of m. 
Using Lemma 16.61 and then Lemma 16.41 we subsequently obtain 


174(5™)! < 


(1 - y/m) 


m 


K (^r7||V,5™r + [||5™r + ||Vp5™r + ||V,5" 

r?||V,5™f+ ^K||Vp5™f 


+ V^K ||V,5™f + ||Vp5^ 


rm||2 


||5™f+ ||Vp5™f+ ||V,5™f +||5^ 


j:m||2 


m ||2 


trm II2 


Finally, choosing r] = r/(m) such that r]{m) —)• 0 and ^ 0, we get that for m 

sufficiently small and for some constant K < oo 


|r4(5™)| = |(C5™,C£^5™)| 


(59) 


< K 


r]{m) + 


r}{m] 


||5^ 


■m 112 I 

m + 


r](m) 


||5^ 
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It remains to consider the cross-term 


T^{6^) = {AC^6^,C6^) + {A6^,CC^5'^) . 
Recalling (|34p we have the following calculations 


= {ACH'^,C5'^) + {A5'^,CC^5^) 

= (1 + y/^)BVp5^ + (1 + - ^h{r){5^ + 1 ) - ^h{r)pVp5^) 

+ (1 + ^)BVr5^ + (1 + y/^)Vrh{r)Vp5^ - ^Vrh{r)p{5'^ + 1 ) 

— \/mh{r)pS/rSi"^) 

= (1 + v^) [{Vp5^,BVr5^) + {\/r6^,B\/p6^)] 


(V, 


, VrbVpS^ 


+ ||v,5" 


-V^[{Vr6^,ph{r)Vpd^) + {Vpd^,ph{r)VrS^)] 

- [{Vp5^, Vrhir)pi6^ + 1)) + {Vr6^, h(r)(5"* + 1))] 
= (1 + V^) [{Vp5^,C^Vr6^) + (V^5™,£^Vp5™)] 

- (1 + v^) [{Vp6^,VpVr6^) + {Vr6^,VpVpS^)] 


(V, 


, VrbVp6^ 




-V^[{Vr6^,ph{r)Vp6^) + {Vp6^,phir)Vr6"^)] 

- [{Vp6^, Vrh{r)p{6^ + 1)) + {VrS'^, h{r){6'^ + 1))] 


Using now Lemma 16.11 on the hrst term of the right hand side of the last display 
we obtain 


r3(j™) = -2(1 + {VpVpS^, VpVrS^) 

- (1 + v^) [{Vp6^, VpVr6^) + {Vr6^, VpVpd”^)] 

{\/pS^,\/rbVp6^) + ||V^<5”^|| 


-V^[{Vrd^,ph{r)Vp5^) + {Vp5^,ph{r)VrS^)] 

- [{Vp6^, Vrhir)pi6^ + 1 )) + {VrS^, hir)i6^ + 1 ))] 

Next we bound terms from above. Using Lemma 16.21 we have for t/ > 0 


173(5™)! < (1 + V^)K ||VpVp5™f + ||VpV,5^ 


1 


1 


T] |Vp5™f + ^ l|VpV^5™f + V l|Vr.5™|r + ^ ||VpVp5' 


Ap 


Ap 


+ 


iL[||V,5™|| ||VpVp5™|| + ||Vp5™|| ||VpV,5™||] 


+ (1 + v^) \k ||Vp5™f + ||V^5™f 
+ [l + ||5™f + ||Vp5™f + ||V^5™f 
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The constant K may change from line to line. Using Lemmas 16.51 and 16.61 we obtain 


173(5™)! < v^(l + ^M)K 1 + ||5™f + ||Vp5™f + ||V^5^ 

rm ||2 


r?|5^ 


+ vm 


+ (1 + y/m)K 
+ (1 + y/m) 


Vfn 

4p 

l + ||5™f+ ||5™||^i] 

m ||2 

\\m 

Ap 

'l + ||5™f+ ||5 


iL||Vp5™f + ||V^5™f 

2 I 11 vT am 112 


+ ^K ||5™|r + ||Vp5™|r + ||V^5" 


Applying then Lemma 16.41 to estimate the term ||Vp5™|p on the fourth line of the 
last display, we obtain the following bound 


1^3(5™)! <(v^ + 


K 


1 + ||5™r + ¥n\m\ + VK WS^fm + (1 + l|V,5^ 


rn\\2 


Finally, choosing r] = r]{m) such that r]{m) —)• 0 and —?> 0, we get that for m 

sufficiently small and for some constant K < oo 


m6n\<{V^+^)K l + ||5™f+ ||5™||^1 +7?(m)A:||5™||^i + (l + V^)||V,5" 
r][m) L J 


This concludes the proof of the lemma. 

Appendix C. Proofs of Lemmas in Section [71 


□ 


Proof of Lemma 7.1 This can be shown by using Theorem 4.2.5 in [T]. Let {Pt)t>o be 
the Markov semigroup corresponding to generator on T- 

By Lemma l6.ll with m = 1, we obtain for the first term (recall that p^{p,r)dpdr 
is the invariant measure corresponding to the operator that the Dirichlet form 
associated with {Pt)t>o can be calculated as follows 

^(/) = f ' f) L'^(y-pO) ~ pf\\L‘^(y-,p°) ■ 

Thus by Theorem 4.2.5 of [T] the validity of Poincare inequality is equivalent to 
exponential convergence to equilibrium of the semigroup {Pt)t>o'- 

j (^Ptf - J {Ptf)p°{p,r)dpdr'^ p^{p,r)dpdr < c{f)e~‘^^/^ . 

for some constant k > 0. The above inequality is true since admits a spectral gap 
(see [To] ). □ 
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Proof of Lemma 1.2. We make use of our equation 
we get 

(60) 


as 


well as Lemma 16.21 and 


1 




- iL||^™3||2,2(y.pO) 


1 


L'^{y\p°) 

for some constant LC > 0 independent of m. 




L^{y-,pO) 


< y/rn\\B^i^2\\L^{y-pO)\\^£^3\\L2(y-pO) , 


have 


We apply Lemma [7.11 using the fact that / "^'^^{p,r)p^{p,r)dpdr = 0, and we 

Jy 


11^ 

< 


m ||2 

i,3\\L^{y;pO) 


\I/"^ _ 


^e,3 ~ 


[ ^e^3iP^r)p°{p,r)dpdT 

Jy 

[ '^e^3iP^r)p°{p,r)dpdT 

Jy 


+ 


LHy;p°) 

2 

LHy;p°) 


ly 


'^t3ip^r)p^{p,r)dpdr 


+ 


+ 


/ '^^3ip,r)p^{p,r)dpdr - / '^^^{p,r)p^{p,r)dpdr 

'y Jy 

2 


^ II 


%\\hiy-,p0) + (^J^'^£^3iP^r)6^{p,r)p^{p,r)dpdi 


< k|| Vr,^y\ll r2('-().„0'i + 11^. 


m ||2 


irm ||2 


£,3\\L‘2{y-,pO)¥''‘ llL2(y;pO) 


' P^£,3\\L^{y-pO) 

Since we have lim ||(^™'||L 2 /y-pO) = 0, we can choose m small enough so that 

m—^O y.J') 

(61) \\^T,3\\L^{y\p°) ^ P^T,3\\L^{y-,p°) ■ 

Combining (j6ip and (j6U|) we see that we have 
1 


(1 — 2KKy/rn) 




L^{y\p°) 


< y/m\\B^l^2\\L'^{y-pO)\\'^T,3\\L'^{y\p° 


) ■ 


Using m again we see that 

2 


(1 — 2KK^/m) 


1 




L^iy-,p0) 


< 2Km\\B'i>ip\\L2^y.pO-) 


Vp^'^3 


L2(y-,p0) 


This means that we have the bound 


1 


Vp'i>T3 


< 


2 Km 


L^{y-,p°) ^ ~ 

Now apply (IHT]) again we obtain the bound 


||^^£,2||L2(y;pO) 




This proves the lemma. 


□ 
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Proof of Lemma 7.3. Let us write 'L in place of for similicity of notations. We set 
/ = and we look for the equation that / satishes: 


C/ 


= lAf + ^Bf 

m \/m 


= —(-P-Vp/ + Ap/) + ^(6(r) -VpZ+p-Vr/) 

m vm 


= —(-P ■ (2Wp^) + 2|Vp^'p + 2^'Ap^') + -^(b{r) ■ 2Wp^' + p • 2^Vr^') 
TTl V TTl 


= 2^(£™^) + —iVp^l 
^ m 


Using the equation (pHjl we see that 


(62) 


CLf = -2V^^B^e,2 + -|Vp^|2 
^ m 


Making use of Lemma 16.11 we have 

{ C ^ fJ ) my , p ^) = -^l|Vp/||i2(y.pO) + ■ P , f ) LHy -, pO ) ■ 

This gives 


(63) l|Vp/||i 2 (y.pO) = -^{h{r) ■ p, f)L^(y.pO) - m{C^f, f)L^(y.pO) • 

Making use of (l6^ . (f63]l . the fact that / > 0 and Lemma f6.2l we get, for some 
constant K > 0 independent of m that may vary from line to line, 


ll^p/llL2(y;pO) 

= ■^^(^(^) ■ P) L‘^{y\p°) + /)L2(y;pO) — 2(|Vp'I'|^, f)l2^y.p0'^ 

< ^{Kr) •P,/^)L2(y;pO) +2m3/2(^-5^-^2,/)L2(y;p0) 

^ ■^^-^ll/llL2(y;pO)||Vp/||j;^2(y.pO) + m^/^(||'I'^'I'£^2|li2(y.pO) + Il/lli2(y.p0)) 

- ■^^-^ll/lli2(y;pO)||Vp/||i2(y.pO) + ?7r^^^(||'I'^|||2(;y.pO)||(^^£,2)^|li2(y.pO) + 11 / 11 L2 (:y ;p0) ) 

(64) 

^ -^[■^^ll/llL2(y;pO)||Vp/||j^2(;y.pO) + 11 / 11 ^2 (y .^0)] . 
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Now we apply Lemma l7. II and we see that for some k > 0 we have 


\\f\\L^(y-,p°) 


< 


/- 

/ f{p,r)p^{p,r)dpdr + / f{p,r 

)p^{p, r)dpdr 



Jy 

Jy 


L'2{y-,p0) 

/- 

/ f{p,r)p^{p,r)dpdr 

+ 

/ f{p,r)p^{p,r)dpdr 


Jy 

L2(y;p0) 

Jy 



< '^l|Vp/llL2(y;pO) + ||'I'||^2(y.pO) ■ 


In the last step we used the fact that / = Now we apply Lemma 17.21 and we 
see that ||'I'||^ 2 (;y.pO) < Krrfi for some constant K > 0 independent of m. Thus we see 
that 

(65) ll/llL2(y;pO) <K[\\Vpf\\l2(^y.pO)+m^] 

Combining (IM)) and (l6^ we see that 

\\'^pf\\l^{y;pO) — -^[V^||Vp/|||2(;y.pO) + 11 Vp/11 ^,2 (y .pO) + 11 Vp/11 ^2 (y .pO) . 

This gives lim ||Vp/||r 2 M; n 0 \ = 0- Apply (165]) again we see that the claim of the 
Lemma follows. □ 


Proof of Lemma \7.f\ The proof of this lemma is completely analogous to that of Lemma 
16.11 and thus it is omitted. □ 
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